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We study the dc Josephson efTect and the density of states in a multiterminal structure of cross- 
type geometry which consists of four superconducting electrodes connected by one-dimensional 
normal or ferromagnetic wires. We find that the Josephson current Ij^ has a sinusoidal dependence 
on the phase difference between the superconductors in the horizontal wire with a critical current 
Ic that can be varied by changing the phase difference <j)y between the superconductors in the 
vertical wire. The period of the function Ijz{(j)z) depends on the ratio of the interface resistances 
Rsn,z/Rsn,y beiug equal to 2-k if this ratio is small and equal to 4-k in the opposite limit. We 
also calculate the amplitudes of both the singlet and the odd-frequency triplet component in the 
system under consideration. The triplet component amplitude may be significantly larger than the 
amplitude of the singlet component. 



I. INTRODUCTION 

In the recent years there has been an increasing interest 
in studies of Josephson junctions (JJs) with controllable 
parameters. The value and sign of the Josephson critical 
current Ic can be changed by variation of temperature 
and external or internal parameters. To some extent, 
this interest is caused by the possibility of technological 
applications of such JJs, for example in Q-bits (for a 
review see e.g. Ref. 1). 

A remarkable example of a JJ with properties which 
depend on internal parameters in a crucial way is 
a superconductor/ferromagnet /superconductor (S/F/S) 
junction. A specific action of the exchange field on the 
spins of Cooper pairs leads, in a certain range of tem- 
peratures and thicknesses of the ferromagnetic films, to 
a sign change of the critical current Ic in S/F/S JJs^. 
This effect predicted long time ago'^'* has been observed 
by many experimental groups^^^^. In S/F/S JJs con- 
sisting of ordinary singlet s-wave superconductors and a 
ferroniagnet with a nonhomogencous magnetization not 
only a singlet component arises in the F due to proximity 
effect, but also a nontrivial, so-called odd triplet compo- 
nent which penetrates the ferromagnetic layer over a long 
distance (for a review see Refs. 15-17). This component 
observed in many experiments^*"^''^ may also lead to an 
unusual and controllable behavior of JJs. 

The latter property has been demonstrated in recent 
experiments^''^ on a multi-layered S/F'FFF'/S Josephson 
junction in which the F' and F layers are a weak (F' 
— >• Ni or Pdo. 88^10.12) and a strong (F-^- Co) ferromag- 
net, respectively. Since the singlet component leaking 
from the superconductor S (S Nb) cannot penetrate 
the rather thick Co layer {dco = 20nm), the Josephson 
coupling is provided by the odd triplet component arising 
due to non-collinearity of magnetization vectors M in the 
F' and F layers. The maximum amplitude of this triplet 
component occurs at perpendicular orientation^^"^''. Ap- 
plying an external magnetic field, the authors of Ref. 24 
changed the mutual orientation of the magnetization vec- 



tors in the F' and F layers and could increase the critical 
current Ic more than by an order of magnitude. 

In S/FF/S JJs where the magnetization vectors M in 
the F layers are parallel or antiparallel with respect to 
each other, the critical current Ic depends on the mu- 
tual orientations of the M vectors (see theoretical papers 
Refs. 28-30). This dependence was measured by Robin- 
son et al."^^ In accordance with the theoretical predictions 
they found that the current Ic is larger for the antiparallel 
orientation than for the parallel one. 

The change of sign of the critical Josephson current 
may be realized not only in S/F/S JJs, but also in 
superconductor/normal metal/superconductor (S/N/S) 
JJs (see theoretical papers Refs. 32-34) if the quasipar- 
ticle distribution function in the normal metal is not in 
equilibrium with the superconducting reservoirs. This 
disequilibrium may be achieved, e.g., in a four terminal 
JJ by passing a dissipative current trough the lateral N 
reservoirs. Then, the critical current Ic{V) decreases as 
a function of the voltage V applied between the N reser- 
voirs and may become negative. This effect has been 
established in experiments by Baselmans et al.'^^ 

The critical current Ic in multi-terminal JJs may not 
only decrease under non-equilibrium conditions, but may 
also increase'^^ . A weak increase of the Josephson current 
has been observed in Josephson weak links under mi- 
crowave radiation several decades ago'^^. The observed 
increase of Ic is probably related to the stimulation of 
superconductivity by microwave irradiation predicted by 
Ehashbcrg'^^. 

Multi-terminal S/N/S JJs were studied in several 
papers'^^''^^''*". These junctions contain a dissipative ele- 
ment, a normal conductor or leads through which a dissi- 
pative current was passed. For example, in experiments 
by Petrashov et al.'^^"^^ the modulation of the conduc- 
tance between the N reservoirs due to a variation of the 
phase difference ip between two superconductors in four- 
terminal structures (such structures are called Andreev 
interferometers) has been measured. 

The possibilities to change the critical Josephson cur- 
rent in non-dissipative four-terminal structures by pass- 
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ing current through lateral electrodes were studied theo- 
retically in Refs. 42 and 43 (ballistic limit) and in a recent 
paper ''^ (dirty case). The authors of Ref. 44 assumed a 
weak proximity effect. Physically, this means that the 
interface resistance is much larger than the resistance of 
the normal two-dimensional region between four super- 
conducting leads. In this limit the Usadel equation can 
be linearized. Solving this linearized equation for a rect- 
angular geometry, Alidoust et al.^"*^ have found the de- 
pendence of the critical current Ic between the left and 
right superconductor on the phase difference between the 
two other electrodes (upper and down). 

In this paper we consider a four-terminal Josephson 
junction of cross- type geometry as shown in Fig. 1. It 
consists of a normal metal (n) wire or film which con- 
nects four superconducting reservoirs. In our analysis 
the connecting wires of the system are considered as 
one-dimensional. The normal metal wires have the total 
lengthes 2Ly [v = z,y) along the corresponding axis v. 
Unlike Ref. 44, we analyze the case of arbitrary strength 
of proximity effect. Using the Usadel equation and as- 
suming that the Thouless energy i^xh = D/L^ is large 
enough [i^Th 3> max(T, A)], we find an exact solution for 
the quasiclassical Green's functions. Here, D = vpl/i is 
the diffusion coefficient, vp the Fermi velocity, / is the 
electronic mean free path, T denotes temperature, A is 
the superconducting gap parameter, L = Lz + Ly, and we 
choose units such that h = ks = 1. Note that the case of 
a high Thouless energy compared to the exchange field 
was analysed in Refs. 45 and 46 for another SF system. 
We consider temperatures T far below the critical tem- 
perature Tc of the superconducting electrodes so that the 
temperature dependence of the gap function A(T) can 
be neglected, i.e. A represents the T = value of A(T). 
With the help of these functions we calculate the density 
of states (DOS) at the crossing point and find the de- 
pendence of the critical Josephson current in the vertical 
JJ on the phase difference in the horizontal JJ. In the 
considered "zero-dimensional" limit expressions for the 
quasiclassical Green's functions are more compact and 
describe the considered JJ for any ratio of the interface 
resistance to the resistance of the n— wires. In particu- 
lar they describe gapped states of the n— wires and the 
dependence of the induced gap on the phase difference. 

We also study the case of ferromagnetic wires connect- 
ing the superconductors (Sec. V). It is assumed that the 
magnetization M in the wires is oriented along the wires 
and that the ferromagnets are not strong so that the con- 
densate penetrates into the ferromagnets and the condi- 
tion i?Tii ^ Vex is fulfilled, where Vc^ denotes the ex- 
change energy. We show that not only the singlet com- 
ponent arises in the ferromagnetic cross, but also the 
odd-frequency s-wave triplet component with non-zero 
projection of the total spin S on the M vectors. The 
amplitude of this triplet component may be significantly 
larger than the amplitude of the singlet component. In 
the concluding section we discuss the obtained results. 




FIG. 1. (Color online) Schematic setup of a four-terminal 
Josephson junction of cross- type geometry consisting of four 
superconducting reservoirs connected by one-dimensional nor- 
mal metal wires. The junction is assumed to be located in 
the zy— plane with interfaces a.t u — C,Lv where v = z,y and 
^ = ±. The superconductors are supposed to be identical but 
with different phases (/p^^ of the order parameter. 

II. MODEL AND BASIC EQUATIONS 

We consider an S/S structure of cross- type geometry as 
shown in Fig. 1. In the horizontal and vertical directions 
we have a S{—Lz)—n/n—S{Lz) and S(— Lj,)— n/n— S(Ly) 
structure, respectively. We will employ the method of 
quasiclassical Green's functions which is well developed 
and can be used in problems where quantum effects on 
the scale of the Fermi wave length are unimportant^^"^^. 
In particular, this theory was applied to the study of 
multiterminal S/S and S/N circuits in many papers (see, 
e.g., Refs. 32, 44, 50-53 and references therein). 

We study the case of dirty superconductors that corre- 
sponds to most experimental situations, i.e. we assume 
that the electronic mean free path / is much shorter than 
the superconducting coherence length = \/ D/ A and 
the lengthes of the n— wires, L^. In the considered equi- 
librium case one needs to find only the retarded (ad- 
vanced) Green's function in order to obtain thermody- 
namical quantities as the density of states or the Joseph- 
son current. For singlet superconductors as assumed here 
these functions have a 2 x 2 matrix structure in the 
particle- hole space. We denote the quasiclassical tem- 
perature Green's function by g{uj) where ui = {2n+ 1)ttT 
is the Matsubara frequency. In the n— wires, the matrix 
function g obeys the Usadel equation^* 

'dAg-d,g) + Kl[T3,9]=0 (1) 
and the normalization condition 

5 -.9 = 1, (2) 

where = ui/D, is the third Pauli matrix, and 1 is 
the 2x2 unity matrix in the particle-hole space. The 
boundary conditions have the form^^ 
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9 ■ d^gl 



Ck^ [g,gs] 



(3) 



where C = ±, — p/2Rns.vi P is the specific resistance 
of the n— wire, RnS,u is the nS interface resistance at 
V = QLi, per unit area, and we assume that RnS.v{Ly) = 

RnS,v(,-Lu). 

As was noted above, the condition 



lI < min{i:>/r,i:>/A} 



(4) 



is supposed to be fulfihed. Therefore, one can integrate 
Eq. (1) over v from — to and from to regarding 
the Green's function g as nearly constant. Taking into 
account the boundary condition, Eq. (3), we obtain the 
equations 
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FIG. 2. (Color online) Density of states A'^(£) in the n— wire 
as a function of normalized energy e/A for = 2n/3 and 
different phase differences 02 . The plots are presented for the 
parameters 7/A = 0.03, = ry = 1, and L/^s = 0.1. 



(g ■ ^i^^L^^o ~ ig^gsiCLi,)] = -^II^ [fg,^] . (5) 

Then, we sum up all these equations and take into ac- 
count the law of " current" conservation 



Jzi+0) + Jy{+0) = M-0) + Jyi-O), 



(6) 



where J^(z) ^ Dg ■ d^g\y^^ and Jy{y) = Dg ■ dyg\^^^. 
We arrive at the equation 



1^1^'^ [^3, .9] + (r^G^ +ryGy] ,g 



= 



(7) 



with 7v = x-uL, Gv = {gs (Lu) + gs (-^i^)] /2. Equation 
(7) can be written in the form 



M,g 



-0, 



(8) 



where M — nf^L'^r^ + {tzGz + fyGy^ . The Green's func- 
tions in the S reservoirs have the form 



M = TaMaN + fs {T2m2 + tiTOi) . 



(11) 



where mi — sin {^z) cos (0z/2)-l-ry sin cos (0^/2), 
m2 = cos (<I)z) cos (02/2) -|- cos ($j,) cos (0j^/2), 
= [kIL'^ + {r^ + Ty) 5s], = + <y3i.-)/2, and 
4'v — Vv+ — H^v-- The solution of this equation is sought 
in the form g = aM^ where the coefficient a can be de- 
termined from the normalization condition, Eq. (2). We 
obtain = [M|^ -|- fgiml -I- to^)] ; where the sign of 
a is chosen in such a way that the density of states is a 
positive quantity. Thus, in the limit of a short Josephson 
junction the solution of the Usadel equation is given by 



M 



M 



(12) 



with Vn = v/M3Ar + /|(™l +^2)- 

Having determined the Green's function g, we can now 
calculate physical observable quantities such as the den- 
sity of states N{e) and the Josephson current Ij. 



gs {CLy) = gsT3 + fs [t2 cos((^,,^) fi sin((pi.^)] , (9) 



III. DENSITY OF STATES 



where Lp^i^ = ip,y{(^L^) is the phase of the order param- 
eter in the superconductor at z = (L^ and y = C^y, 
respectively. For simplicity we assume that the super- 
conductors are made of identical materials, i.e. their su- 
perconducting gap parameters A have equal magnitudes 
and therefore the functions 



gs 



fs^^, = Voj' + A2 



(10) 



are identical for all superconductors. Using Eqs. (9) and 
(10), the matrix M can be represented in the form 



With help of the solution for the Green's function g, 
Eq. (12), we can find the density of states N{e) in the 
n— wire, 



N{e) = Re 



(13) 



One can consider the following limiting cases of Eq. (13). 

For very high interface resistances RnS,^ ^ Rl (i-e. 
Tu ^ 1), where i?^ = Lp/2, one can neglect the term 
containing fs in the definition of so that A^(e) « 1. 
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For very low interface resistances RnS.v ^ 
(i.e. Ti, 3> 1) we set for simplicity 
Ky ~ K, i.e. fz — ry ^ r. Then, M^n 



Rl 



2rgs and Vn « 2r[gl 



ifs/2f{f^l + ^,l)Y/^ with 
/ii = sin($2)cos(02/2) + sin($^y)cos(0;,/2), ^2 = 
cos($2) cos((^z/2) + cos($y) cos(0y/2). This results in 

we obtain a 



N{e) « Re 



I.e. 



/£2_A2 (M?+Mi)/4 
gap induced in the n— wire equal to A in the absence of 
superconducting phases ($y = 0^ = 0) and turning to 
zero for /ii = /i2 = 0. 

In the limit of small energies (e/A ^ 1) 
and small values of r we obtain N{e) w 

Re ' ' 



At 



0, 



2, this expression coincides with 



i.e. at /^i = and /i2 = 
McMillan's formula^^. 

In Fig. 2 we show the dependence of the density of 
states N{e) on the normalized energy e/A for different 
phase differences (f>z keeping (j)y fixed at 27r/3. Here, 
we have set (p„+ = —(pv- ~ (j)u/2 so that $y = and 
TOi = 0. Moreover, since the DOS is calculated from 
N{e) ~ Re [fa • g^] , i.e. from the retarded Green's func- 
tion, wc replace in Eq. (13) e/A ^> (e + i7)/A and choose 
a damping 7/A = 0.03 according to analyticity of in 
the upper half-plane. 

As already mentioned above we find a gap induced in 
the n— wire. The position of the gap with respect to 
normalized energy e/A can be tuned by variation of the 
parameters and the superconducting phase differences 
(j)^. In particular, it can be seen that for fixed phase 
difference 4>v one can always find a phase difference = 
2 arccos[— T-p cos(0p /2) /r^] such that the density of states 
becomes unity for all energies, i.e. the gap turns to zero. 
Here, D = y,z denotes the axis perpendicular to = z, y, 
i.e. if for instance v = z then v = y. Since A^(e) is 
proportional to cos((/)2/2) due to symmetry it is sufficient 
to consider phase differences 4>z G [0, 27r] for constant 0^. 

The dependence of the DOS on the phase difference be- 
tween the superconducting leads in S/N/S systems has 
been studied in many papers (see, e.g., Ref. 57 as well 
as Ref. 58 and references therein). Indeed, the proxim- 
ity induced gap in the n— wire has also been obtained 
in Ref. 58 for a three-terminal diffusive S/N/S junction 
with normal metal wires. In agreement with our results, 
the authors of Ref. 58 report a proximity induced gap in 
the DOS of the n— wire which depends on the supercon- 
ducting phase difference (/>. 
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FIG. 3. (Color online) Normalized supercurrent j{Lz) in the 
horizontal wire at the right interface as a function of the hori- 
zontal phase difference 02 for different values of the resistance 
ratio ry = Rl / RnS,y The plots are presented for the two 
cases (a) if^+ = —ipv- = 0i,/2 [^i, = 0] with v = z,y and 

(b) = -if^- = (t)^/2 = 0], •Py+ = = ^y [(jiy = O] . 

The parameters are = 1, T/A = 0.05, L/^s = 0.1, (a) 
<j)y = 27r/5, and (b) $y = 27r/5. 



where Tr is the trace over particle-hole indices. Sub- 
stituting the expressions for g and gs, Eqs. (12) and 
(9), into this expression, we obtain for the normalized 
supcrcurrents j{CLu) = IjiC^u) /Ijou [where Ijq^ = 
7rA/(ei?„s,i.)] 



T 



3{C,Ly) "^-^[-j sin {4>^) + Cr-p cos ((?!)p/2) 

A2 



X sin {if^c - X! 



1 



P + A2 Vn 



(15) 



IV. JOSEPHSON CURRENT 

By means of Eq. (12) we can also calculate the Joseph- 
son current flowing through the four terminals located at 
V = (L^ given by 

IjiCL.)^Cj^^J2^r{r,-[g,gs]} , (14) 



This formula determines the Joscphson current between 
identical superconductors (a generalization to the case 
of different superconductors is straightforward) for ar- 
bitrary phases of the order parameter and arbitrary 
strength of proximity effect, i.e. the parameters are 
not restricted to particular limits. It is obvious that, as it 
should be, a shift of all the phases by a constant does not 
change the currents j{CL^). In the following we consider 
two particular cases with respect to the superconducting 
phases Lp^^^. 
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FIG. 4. (Color online) Normalized supercurrent j{Lz) in the horizontal wire at the right interface as a function of the 
horizontal phase difference (f}^ for different values of the vertical phase difference (f)y. The plots are presented for the case 
ipv+ = —ifv- = 4'v/'i [$1/ = 0] with u = z,y and the parameters T/A = 0.05, L/^s = 0.1, (a) = 1, = 10, (b) = ry = 1, 
and (c) = 10, ry — 1. 



First, we assume that ip„+ = —(p^- = so that 

$1/ = 0. Then, one can easily obtain from Eq. (15) 
the Josephson current in the horizontal and vertical 
branches. 



-sm( 
X sin ( 



+ cos ((^p/2) 

A2 1 



a;2 + A^Vn' 



(16) 



One can see that if RnS,u ^ RnS,u, i-e. r^, ^ rp, the 
dependence of the Josephson current i{L^) on the phase 
difference 4>i, is close to the sinusoidal one with the period 
2tt (note that in the considered particular case mi ~ 0, 
but the function TO2 ^ and it also depends on 0jy). 

In the opposite limit RnS.v ^ RnS.u, i-e. Ti, <^ rp, 
the phase dependence of the current j{L^) is changed 
drastically as it is almost sinusoidal but the period of 
oscillations is equal to Air. 

The obtained phase dependence of the Josephson cur- 
rent leads to an interesting phenomenon. As is well 
known (see, for example, Refs. 59-61), if a Joseph- 
son junction is irradiated by microwaves with frequency 
Lo and if a bias voltage V is applied between the su- 
perconducting leads the so-called Shapiro steps appear 
on the current- voltage characteristics I{V). The posi- 
tions of the Shapiro steps are determined by the relation 
Vsh.n ~ huj /2en with n = 1, 2, . . . Thus, the main Shapiro 
step corresponds to the voltage Vsii,i = fiaj/2e. However, 
given a bias voltage V is applied between the supercon- 
ducting leads in the horizontal wire and the phases are 
varying in time with frequencies much less than the char- 
acteristic energies (A or T^) of the problem, one can con- 
clude from Eq. (16) that for <^ ry the positions of 
the Shapiro steps are changed: Vsh,n(?'2 ^ r^) = hw/en. 
Moreover, the positions of the Shapiro steps depend on 
the phase difference (j)y, that is, on the current flowing 
through the vertical wire. 



Note that the Att periodic phase dependence of the 
Josephson current is typical for junctions with Majorana 
bound states^^"^^. 

Second, consider now the case when in the horizontal 
branch tpz+ = —fz- = </>z/2, i.e. there is no change 
compared to the previous case, but in the vertical branch 



jiCLz) 



0). Then, we find 



^ sin ((/),) + Ty sin (0,/2 - Q^y) 



A L 2 

xE 



A2 



(17) 



and 



j{CLy) = C-r, cos (0,/2) sin {^y) 
A2 1 



'A 

51 ■,,,2 ^ A2 



w2 + A2 Vn 



(18) 



In this case the currents through the interfaces at z ~ 
C,Lz are different if Vy ^ r^. The currents through the 
interfaces at y = C,Ly flow in opposite directions. 

In Fig. 3 we plot the dependence of the current j{Lz) 
on the horizontal phase difference (/>2 for different values 
of the vertical resistance ratio ry and fixed horizontal 
resistance ratio r^ ^ I. The curves are presented for the 
cases ip,y+ = —(pu- = 4>v/'^ with 4>y = 27r/5, Fig. 3(a), 
and Lpz+ = -Pz- = 0^/2, ipy+ = (py- = $y = 27r/5, 
Fig. 3(b). 

From Fig. 3(a) it can be nicely seen how the phase 
dependence of the supercurrent j{Lz) changes from a An 
periodicity for ry r^ (dotted curve) to a 27r periodicity 
for ry <g; rz (solid curve). In the intermediate region 
where the parameters r^, and ry are of the same order 
of magnitude (dashed curve) we find two local extremes 
close to (bz = 27r. 
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FIG. 5. (Color online) Normalized critical supercurrent 
jc{Lz) in the horizontal wire at the right interface as a func- 
tion of the phase difference 0^ for different values of the re- 
sistance ratios where v = z,y. The plots are presented 
for the two cases (a) (pv+ = —^fiv- = (fiv/i [^v = 0] and (b) 

Vz+ = = 0z/2[<l>2 = 0], iPy+ = iPy- = [(?!>y = O] . 

The parameters are T/A = 0.05 and L/^s = 0.1. 



Prom Fig. 3(b) we conclude that for ^ (solid 
curve) there also emerges a 2tt periodicity that, compared 
to Fig. 3(a), is slightly shifted in the horizontal direction 
due to the different dependencies of the second term in 
Eqs. (16) and (17) on the phase difference 0j, and the 
phase sum $y, respectively. Decreasing the parameter 
ry further results in a perfect sinusoidal dependence on 
(j)z- Here, too, in the limit ^ Vy (dotted curve) the 
periodicity of the supercurrent j{Lz) is A-k but as for the 
case Tz 3> Ty it is shifted in the horizontal direction. For 
the intermediate region ~ Vy (dashed curve) we find 
a different influence of the resistance ratios Ti, than in 
Fig. 3(a). In this case we find a double peak rather than 
two local extremes. 

In Fig. 4 we show the dependence of the supercurrent 
j{Lz) in the horizontal wire at the right interface on the 
horizontal phase difference 4)z for different values of the 
vertical phase difference (j)y and the resistance ratios r^. 
Here, we only analyze the case (p,^+ = —ipv- = 4>v/'2- 

As can be seen from Fig. 4(a), where the case <C Vy 
is displayed, the supercurrent j{Lz) has a perfect sinu- 
soidal shape with a 47r periodicity for arbitrary vertical 
phase differences (j)y. This can be understood by means 
of Eq. (16) where for <^ ry the first term in the square 



ry and four 
In contrast to 



brackets can be neglected compared to the second term 
so that the current is proportional to cos(0y/2) sin(02/2). 
Therefore, due to the dependence on cos(<^y/2), the vari- 
ation of the vertical phase difference (f>y from to 2tt 
leads to a sign change oi j{Lz). 

In the opposite limit r^ S> ry, which is shown in 
Fig. 4(c), the first term in the square brackets of Eq. (16) 
dominates with respect to the second term, i.e. in this 
case the supercurrent j{Lz) is proportional to sm{(l)z) and 
is independent of the vertical phase difference (py. Thus, 
the curves collapse into a single curve with periodicity 
2tt. 

In the intermediate regime w ry the form of super- 
current j{Lz) can be considerably changed by variation 
of the vertical phase difference (jiy . While the dependence 
of the current on the horizontal phase difference (pz for 
0y = 0, 27r is nearly sinusoidal, we find a clearly distinct 
behavior for <f>y = 27r/3,47r/3 where the dependence on 
<f>z is more complicated. 

Finally, we calculate numerically the normalized criti- 
cal Josephson current jc{Lz) in the horizontal wire as a 
function of the horizontal phase difference (pz for different 
values of the resistance ratios r^, . For each set of parame- 
ters r^, we calculate for each phase difference (j)z € [0, 47r] 
the maximal supercurrent ic{Lz) as a function of (fiy and 
plot it versus 4>z as shown in Fig. 5. Interestingly, we 
obtain for the case ipy+ = —(p^- = '/'i//2, Fig. 5(a), two 
maximal values of jc{Lz) for the case r 
maxima of jc{Lz) for the case 3> ry 
Fig. 5(a), in the case where we have in the vertical wire 
ipyj^ = Lpy^ = $y. Fig. 5(b), wc always find four maxi- 
mal values of jc{Lz) and there does not exist a horizontal 
phase difference (f>z for which the critical current vanishes. 

In Fig. 5 the number of peaks of the critical current 
reflects the periodicity of the Josephson current with 
respect to (j)z which can be understood by considering 
Eqs. (16) and (17), respectively. Considering the up- 
per panel of Fig. 5 we find for the case ^ ry (solid 
curve) two peaks of the critical current because accord- 
ing to Eq. (16) the prefactor of the Josephson current 
j{Lz) in this limit is proportional to sin(02/2), i.e. it 
has a 47r periodicity [also see Fig. 4(a)]. On the contrary, 
for rz » ry (dotted curve) we obtain j{Lz) oc sm{(t)z), 
i.e. the Josephson current has a 27r periodicity result- 
ing in four peaks of the according critical current. Sim- 
ilarly, we can explain the appearance of four peaks in 
Fig. 5(b). Here, in the limit ^ ry (solid curve) the 
Josephson current is proportional to sm{(j)z/2 — $-y) and 
in the opposite limit (dotted curve) it is proportional 
to sui((j)z). Therefore, in the latter case we expect four 
peaks of jc{Lz) whereas for ^ ry the critical current 
approaches a constant as a function of (jjz because, here, 
the superconducting phase difference (j)^ only leads to a 
horizontal shift of the Josephson current and docs not 
affect its magnitude. 

According to Eq. (4), where the condition <C 
D/A = ^1 is supposed to be fulfilled. Figs. 3-5 are pre- 
sented for a small value of normalized junction length 
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L/£,s which is kept fixed at = 0.1. 



V. FERROMAGNETIC WIRES WITH 
NON-COLLINEAR MAGNETIZATION 

In this section we consider the case when the cross is 
formed not by normal wires but by ferromagnetic wires 
witli non-collinear magnetization vectors M. We assume 
that the vector M is ahgned in the z— direction in the 
horizontal ferromagnetic wire and along the y— direction 
in the vertical ferromagnetic wire, i.e. in both cases M 
points along the corresponding wire. In order to describe 
the condensate in this case the quasiclassical Green's 
functions should be generalized. They become 4x4 ma- 
trices g in the product space of particle-hole and spin 
variables (see, e.g., Ref. 15). Equation (1) can be writ- 
ten as 



~dz{g ■ dzg) + k^[t3 • aa,g] - «K|^[f3 ■ 0-3, 5] = (19) 
-dy{g ■ dyCj) + K^[f3 • (7o,g] - iK|[fo • 0-2,5] = 0, (20) 

where n^p = Ve^/D and Vox is the exchange energy. The 
matrices fi and ai operate in the particle-hole and spin 
space, respectively (Pauli matrices for i = 1,2,3 and 
identity matrices for i = 0). We assume that the bound- 
ary condition, Eq. (3), is valid in the generalized form 
(5 ~^ 9t gs ^ 9s) also in the considered case, i.e. we 
assume that the S/F interfaces do not affect the spins. 

Following the procedure presented in Sec. II, we obtain 
a generalized version of Eq. (8), 



One can exclude the term proportional to fo • (T2 in 
the matrix Mp by making a transformation (see Ref. 15) 
that describes a rotation in the spin space. 



Ml 







U'^ -g- U, (25) 



where the transformation matrix depends on a via 
C7(a) = exp(if3 • aia/2) = cos(a/2) -I- isin(a/2)f3 • cti. 
One can see that the matrices Gu remain unchanged 
(they commute with U) because they describe the sin- 
glet superconductors which serve as reservoirs and are 
not affected by a rotation in the spin space. Similarly, 
the first term in Mp^a) commutes with U . Therefore, 
after this transformation Eq. (21) attains the form 



[A^F,0] =0, 
where M.p ^ Mp{Q) is given by 



(26) 



Mf = ntL^ T3 • CTo - WcxT3 • o's + rzGz + TyGy. (27) 

Due to the diagonal structure ol M.p in the spin space 
the equations for different spin sectors decouple, and the 
diagonal elements Qs {s = ±, 0+ = gn, g_ = §22) of g 
can be obtained as it was done in Sec. II. This yields 



M 



Fs 



Vr. 



(28) 



with 



[MF,g\ =0 



(21) 



with Mp given by 

Mp{a) = K^L^ T3 • (To - ivox [sin(Q;)'ro • (T2 4- cos(a)f3 ■ (T3] 



-ryGy 



TzGz 



(22) 



where Vc 



K^pL^I Li + L"^ is the normalized ex- 



change energy, tan(a) = 
[gs {Lu) + gs {-Lu)] /2 with 



Ly/Lz, and 



Mfs 
Vfs 



T3 (M^N - siuex) + sfs (f2m2 -f fimi) (29) 



(MsN - sivcxf + fg {ml + ml) 



(30) 



satisfying 2?i?- = T^*f+ for real u. The resulting matrix 
Q can be written as 



00 • CTQ + 03 • 0-3, 



(31) 



where §0,3 = (0+ ±0-) /2. Exploiting the representation 
oi M.F in terms of its diagonal blocks jiA ps as 



gs{C,Ly) = gsh ■ CTo 

+ fs [t2 cos{ip^q) + Ti sm{ip^()] ■ (73. (23) 

As before, Eqs. (19) and (20) must be solved together 
with the normalization condition 

g-g^i. (24) 

At Vex = the matrix element (5)11 in the spin space 
coincides with the matrix g given by Eq. (12). 



Mf^Y.-^p-- ^^^0 + ''^^^ (32) 

S 

WG arrive at 

q = Mf (Rc [Dpi] ao + i Im [Dpi] a^) . (33) 

Using Eq. (25) one can obtain the matrix g = U ■ g ■ W . 
Since A4 f simply transforms back to Mp it is sufficient 
to transform (T3 and we obtain 



8 



where Tr is the trace over particle-hole and spin indices. 



g = Mf(Rc[P^y ao (34) 
+zlm [Dp^.] {cos(a)fo ■ 0-3 + sin(Q;)f3 • ^2} 

Note that the real and imaginary parts of 2?^^ present 
in the equations above must be evaluated treating uj as 
real before applying analytical continuation via, e.g., 



We are interested in the condensate function / (the 
part of g non-diagonal in the particle-hole space) which 
has the form f = f^g + /tr with the singlet component 



fsg = Re [Dp^] fs (f2m2 -f- fimi) ■ 0-3 



(36) 



N{e) = Re 



iTr{(f3 -CTo) ■ g} 



(35) and the triplet component 
I 



/tr = i Im [Dp]] fs {[t2 ■ a-Q cos(a) + ti ■ &i sm{a)] ms + [fi • &o cos(a) - fa • (Ti sin(Q;)] mi} 



(37) 



In Eq. (37) the terms proportional to (Tq describe the 
triplet component with zero projection of the total spin 
of the Cooper pair on the z— axis (/tr,o) whereas the 
terms proportional to ai describe the component with 
non-zero projection of the total spin of the Cooper pair 
on the z— axis (/tr,i)- The appearance of this "Sz 0"- 
triplet component which corresponds to a long-range 
triplet component in long ferromagnets (see Ref. 15) is 
easily understood from the physical point of view. In the 
absence of the horizontal wire, only singlet and triplet 
components with zero projection on the y— axis {Sy = 0) 
are induced in the vertical wire due to proximity effect. 
On the other hand, this triplet component has a nonzero 
projection on the z— axis (5*2 ^ 0). 

In Fig. 6 we plot the ratio of the " Sz 7^ 0" -triplet com- 
ponent amplitude to that of the singlet component, i.e. 
6tr/sg I sin(Q;)Im[P^^]/Re[Pj^^]|, as a function of nor- 
malized energy e/A for the case = ^</'iy+ = 0iy/2 
(<f>^ = 0) and for the normalized exchange energies 
Kx/A 5 [Fig. 6(a)], 1/cx/A = 20 [Fig. 6(b)]. Moreover, 
the figure displays the ratio ^ti/sg for different phase dif- 
ferences (j)z and fixed phase difference (j)y = 27r/3. As can 
be seen from Fig. 6(a), for small normalized exchange 
energies V^:^/ A the ratio 6tr/sg is peaked where the peak 
position depends on the values of the superconducting 
phase differences and the resistance ratios r^. Here, 
the ratio 6tr/sg is smaller than one for the whole energy 
range e/A < 1. Interestingly, the situation changes dras- 
tically upon increasing the normalized exchange energy 
Vex/ A as shown in Fig. 6(b). In this case, the ratio 6tr/sg 
may be significantly larger than one for small normalized 
energies e/A. This can be easily understood by consid- 
ering in Eq. (30) the limit IM3JV — iwoxl ^ /s v^^if + mf , 
which corresponds to the limit l^x/A ^ 1, so that 
btr/sg = I sin(a)ucx/A^3Ar| oc I sin(a)K!x/e|- Thus, we ob- 
tain in the limit e/ A — > for large normalized exchange 
energies a ratio 6tr/sg that may be larger than one. 

In Fig. 6(b), where we have chosen Vex/ A ~ 20, the 
ratio 6tr/sg turns to zero for small normalized energies 



e/A and the phase differences 0z = (solid line) and 
4>z — 27r/3 (dash-dotted line) because in this case the 
normalized exchange energy T4x/A is still too small to 
shift the ratio b^i/sg in this energy range to values signif- 
icantly larger than one. Effectively, we observe that the 
curves in Fig. 6(b) approximately converge to a single 
one with the energy dependence 6tr/sg | sin(a)14x/e| 
for really large normalized exchange energies V^x/A ^ 1. 
Moreover, for large enough values of X^x/A the function 
^tr/sg is always larger than one for the relevant normal- 
ized energies e/A < 1, i.e. for large normalized exchange 
energies the triplet component always dominates the sin- 
glet component. 

Note that the DOS and the Josephson cur- 
rent are determined by the f3(To-component 
530 = (Re[2?^^]M3Ar -I- Im[X»^]|^]we3,)f3 • (To and the 
singlet component fsg of g, respectively. In particular, 
the Josephson current is determined by Eq. (15) in 
which the factor VJ^^ should be replaced by Re[r']^]J_]. 

It is important to emphasize a difference between a 
long and short F-wire in the presence of the singlet and 
triplet components. As is known-'^^'^^, th e single t compo- 
nent decays over a short length ^p ^ y/DjVx from the 
S/F interface, whereas the triplet component penetrates 
into the F wire over a much longer length which does not 
depend on the exchange energy Vex- Near the S/F inter- 
face the charge is carried out by the singlet component 
because the commutator at the right-hand side of Eq. (3) 
turns to zero due to "orthogonality" of the singlet and 
triplet components. At distances |z| ^ S,p with respect to 
the interfaces, a conversion of singlet and triplet currents 
occurs and at distances |z| ^ £,p with respect to the in- 
terfaces the current is due to the triplet component. In 
short S/F/S junction such a conversion does not occur 
and the current is only caused by the singlet component. 

In Fig. 7 we show the DOS in the n— wire as a function 
of normalized energy e/A for the case ipi,^ ~ —^u- = 
(\>y/2 ($1^ = 0) with 0z = 27r, (\>y — 27r/3 and for different 
exchange energies Vex/ A. 
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<Pz 


= 


-- (t>z 


= tt/2 


■■■■ <l>z 


= 3tt/2 


■-■ <j>z 


= 2ti 



<t>z 


= 


-- (t>z 


= 7r/2 


■■■■ <Pz 


= 37r/2 


■-■ <j>z 


= 271 



"0 0.2 0.4 0.6 0.8 1 

e/A 

FIG. 6. (Color online) Ratio of the "5^ / 0" -triplet compo- 
nent amplitude to that of the singlet component as a function 
of normalized energy e/A for the case ^p^^ = —ipv- = 0,^/2 
= 0), (j)y = 27r/3, and different phase differences (p^. The 
plots display the function btr/sg for 7/A — 0.03, = Vy = 1, 
L.lis = Lj,/Cs = 0.1, and (a) Kx/A = 5, (b) Kx/A = 20. 



We observe that the effect of the normalized exchange 
energy Vi^x/A on the DOS is a splitting of the peak for 
vanishing T4x/A, i.e. the gap induced in the n— wire is 
shifted to lower normalized energies. Increasing the nor- 
malized exchange energy Vox/A the two peaks initially 
move further apart from each other and their amplitudes 
decrease. If the peak located at lower normalized energy 
approaches the origin by increasing T4x/ A furthermore, 
it bounces back from the vertical axis and also moves to 
higher energies together with the other peak. Thus, once 
the peak at lower normalized energy has bounced back 
from the vertical axis, the induced gap in the n— wire van- 
ishes and the DOS approaches 1/2 between the two peaks 
and 1 at the tails. Finally, as the normalized exchange 
energy is increased further the region where the DOS ap- 
proaches 1/2 becomes narrower and the two peaks move 
to higher energies. 

This feature has also been obtained in Ref. 66 where a 
diffusive ferromagnet attached to a normal metal and a 
superconducting reservoir was considered. It is reported 
that for certain values of exchange field a peak at zero en- 
ergy appears in the DOS which corresponds to the gener- 
ation of the triplet pairing amplitude at low energy. The 
authors elaborated that the form of the DOS in the fer- 




Kx/A = 

Vex /A = 10 
Vex /A = 20 
^ex / A = 40 



FIG. 7. (Color online) Density of states A^(£) in the n— wire 
as a function of normalized energy e/A for the case i^uj^ — 
— tfiv- = 4>v/2 ($1, = 0), (f)z = Stt, 4>v = 27r/3, and different 
exchange energies Vex/ A. The plots are presented for the 
parameters 7/A = 0.03, — Vy = 1, and L^/^^s = Ly/^s = 
0.1. 



romagnet for low energies can be used to find the pairing 
symmetry in the ferromagnet, i.e. whether the peak is 
generated by the singlet or triplet pairing states. From 
this analysis we can conclude that since in Fig. 6(b) for 
4>z = 27r (dash-dotted line) ) the triplet component dom- 
inates the singlet component the peak at zero energy in 
Fig. 7 for 14x/A = 20 (dotted line) is generated by the 
triplet component at low energies. 

Note that the influence of the triplet component on 
the DOS was also analyzed in Refs. 46, 67, and 68 for 
other systems. Also, the splitting of the main peak in 
the DOS due to the exchange field was studied earlier in 
Rcfs. 69-72. 

In order to figure out the influence of the normal- 
ized exchange energy 14x/A on the supercurrent, in 
Fig. 8 we plot the normalized Josephson current at the 
right horizontal interface, i.e. at z = iz, for the case 
tf^+ = — = 01//2 ($,y = 0) as a function of the 
phase difference (f>z with the parameters i^j, = 27r/5 
Lz/^s = Ly/^s = 0.1, - 1, = 0.1 [Fig. 8(a)], r„ = 1 
[Fig. 8(b)], and for different values of the exchange energy 
14x/A. 

Here, we find that the effect of the normalized ex- 
change energy T4x/A on the supercurrent is a distor- 
tion of the sinusoidal shape while the 27r periodicity for 
rz = I, ry = 0.1 is conserved [Fig. 8(a)]. With increas- 
ing normalized exchange energy the supercurrent is more 
and more distorted leftwards together with a decrease in 
amplitude, i.e. the exchange energy can significantly re- 
duce the magnitude of the supercurrent. For the case 
Tz = fy = 1, Fig. 8(b), the effect of the normalized ex- 
change energy is weaker compared to the case = 1, 
Vy = 0.1, Fig. 8(a). Here, much higher values of Kx/A 
are necessary in order to reduce the magnitude of j{Lz) 
considerably but the local extremes close to (pz = 2tt are 
smoothened more easily. For the case <^ Vy the form 



of the supercurrent is nearly unchanged by the normal- 
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Vex/A = 
Vex/A = 20 
Vex/A = 40 



— Ve,/A = 
--Vex/A = 20 
■■■■Vex/A = 40 



FIG. 8. (Color online) Normalized supercurrent j{Lz) in the 
horizontal wire at the right interface as a function of the hor- 
izontal phase difference (j>z for different values of the nor- 
malized exchange energy Kx/A. The plots are presented 
for the case ip^+ = — (/?^- = <j},^/2 ("1?^ = 0) with v = z, y 
and the parameters T/A = 0.05, = 1, 0y = 27r/5, 
L,/is = Ly/is = 0.1, (a) = 0.1, and (b) Vy = 1. 



ized exchange energy V^cx/A and retains its sinusoidal 
behavior with a Air periodicity [see Fig. 3(a)]. 

Comparing the critical Josephson current jc{Lz) at 
the right horizontal interface with respect to the normal 
metal wire and the ferromagnetic wire, we find that the 
magnitude of jdLz) decreases with increasing normal- 



ized exchange energy T4x/A while the form of the critical 
Josephson current remains unchanged (see Fig. 5). Here, 
the decreasing effect of the exchange energy is barely vis- 



ible for Vz 
for r-, = r 



;S> Ty and Tz 



<C Ty while it is slightly stronger 



VI. CONCLUSION 

We have calculated the Josephson dc current in a mul- 
titerminal structure of cross-type geometry which con- 
sists of four superconductors and one-dimensional nor- 
mal or ferromagnetic wires connecting the superconduct- 
ing electrodes. The length of the wires is assumed to be 
short, i.e. the Thouless energy exceeds any characteristic 
energy of the system (A, T, or 14x). In this case one can 
easily obtain a solution of the Usadel equation for the 
quasiclassical matrix Green's functions ^ in a compact 
form. With the help of the matrix g we find the DOS 
and the dc Josephson current Ij. It turns out that both 
the DOS and Ij strongly depend on the phase differ- 
ences between opposite superconductors. The sinusoidal 
dependence of Ij on the phase difference 0^ has different 
period which is determined by the ratio of the interface 
resistances Rsn,z/Rsn,y If this ratio is small, the period 
is equal to 27r, in the opposite limit the period equals An. 

It is found that the DOS in the wire has a minigap 
which depends on the phase difference and can turn 
to zero. 

We have obtained also simple formulas for the sin- 
glet and triplet components induced in the ferromagnetic 
wires. The odd-frequency s-wave triplet component has 
non-zero projections of the total spin S on the directions 
on the wires and can be larger than the amplitude of the 
singlet component. However, the Josephson current is 
caused only by the singlet component. 

The obtained results may be useful in superconducting 
electronics and spintronics (see Ref. 44 and references 
therein). 



^ Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 

73, 357 (2001). 
^ A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 
^ L. N. Bulaevskii, V. V. Kuzii, and A. A. Sobyanin, JETP 

Lett. 25, 290 (1977). 
* A. I. Buzdin, L. N. Bulaevskii, and S. V. Panyukov, JETP 

Lett. 35, 178 (1982). 
^ V. V. Ryazanov, V. A. Oboznov, A. Y. Rusanov, A. V. 

Veretennikov, A. A. Golubov, and J. Aarts, Phys. Rev. 

Lett. 86, 2427 (2001). 
« V. A. Oboznov, V. V. Bol'ginov, A. K. Feofanov, V. V. 

Ryazanov, and A. I. Buzdin, Phys. Rev. Lett. 96, 197003 

(2006). 

T. Kontos, M. Aprili, J. Lesueur, F. Genet, B. Stephanidis, 
and R. Boursier, Phys. Rev. Lett. 89, 137007 (2002). 



* Y. Blum, A. Tsukernik, M. Karpovski, and A. Palevski, 

Phys. Rev. Lett. 89, 187004 (2002). 
^ A. Bauer, J. Bentner, M. Aprih, M. L. DeUa Rocca, 

M. Reinwald, W. Wegscheider, and C. Strunk, Phys. Rev. 

Lett. 92, 217001 (2004). 
" H. Sellier, C. Baraduc, F. Lefioch, and R. Calemczuk, 

Phys. Rev. Lett. 92, 257005 (2004). 
^ J. W. A. Robinson, S. Piano, G. BurneU, C. Bell, and 

M. G. Blamire, Phys. Rev. Lett. 97, 177003 (2006). 
^ M. Weides, M. Kemmler, H. Kohlstedt, R. Waser, 

D. Koelle, R. Kleiner, and E. Goldobin, Phys. Rev. Lett. 

97, 247001 (2006). 

A. A. Bannykh, J. Pfeiffer, V. S. Stolyarov, L E. Batov, 
V. V. Ryazanov, and M. Weides, Phys. Rev. B 79, 054501 
(2009). 



11 



D. Spmngmann, K. Westerholt, H. Zabel, M. Weides, and 

H. Kohlstedt, J. Phys. D 42, 075005 (2009). 

F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Rev. Mod. 
Phys. 77, 1321 (2005). 

1® M. Eschrig, Physics Today 64, 43 (2011). 

K. B. Efetov, I. A. Garifulhn, A. F. Volkov, and K. West- 
erhoh, in Magnetic Nano structures: Spin Dynamics and 
Spin Transport, edited by H. Zabel and M. Farle (Springer, 
Berlin, 2012). 

R. S. Keizer, S. T. B. Goennenwein, T. M. Klapwijk, 

G. Miao, G. Xiao, and A. Gupta, Nature 439, 825 (2006). 

I. Sosnin, H. Cho, V. T. Petrashov, and A. F. Volkov, Phys. 
Rev. Lett. 96, 157002 (2006). 

2° T. S. Khaire, M. A. Khasawneh, W. P. Pratt, and N. O. 
Birge, Phys. Rev. Lett. 104, 137002 (2010). 
J. W. A. Robinson, J. D. S. Witt, and M. G. Blamire, 
Science 329, 59 (2010). 

D. Sprungmann, K. Westerholt, H. Zabel, M. Weides, and 

H. Kohlstedt, Phys. Rev. B 82, 060505 (2010). 

M. S. Anwar, F. Czeschka, M. Hesselberth, M. Porcu, and 

J. Aarts, Phys. Rev. B 82, 100501 (2010). 

C. Klose, T. S. Khaire, Y. Wang, W. P. Pratt, N. O. Birge, 

B. J. McMorran, T. P. Ginley, J. A. Borchers, B. J. Kirby, 

B. B. Maranville, et al., Phys. Rev. Lett. 108, 127002 

(2012). 

'^^ M. Houzet and A. I. Buzdin, Phys. Rev. B 76, 060504 
(2007). 

2^ A. F. Volkov and K. B. Efetov, Phys. Rev. B 81, 144522 
(2010). 

L. Trifunovic and Z. Radovic, Phys. Rev. B 82, 020505 
(2010). 

F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. 
Lett. 86, 3140 (2001). 

Y. M. Blanter and F. W. J. Hekking, Phys. Rev. B 69, 
024525 (2004). 
3" A. Zaitsev, JETP Lett. 90, 475 (2009). 

J. W. A. Robinson, G. B. Halasz, A. I. Buzdin, and M. G. 
Blamire, Phys. Rev. Lett. 104, 207001 (2010). 
A. F. Volkov, Phys. Rev. Lett. 74, 4730 (1995). 
S.-K. Yip, Phys. Rev. B 58, 5803 (1998). 

F. K. Wilhelm, G. Schon, and A. D. Zaikin, Phys. Rev. 
Lett. 81, 1682 (1998). 

J. J. A. Baselmans, A. F. Morpurgo, B. J. van Wees, and 
T. Klapwijk, Nature 397, 43 (1999). 
'•^^ R. Seviour and A. F. Volkov, Phys. Rev. B 61, R9273 
(2000). 

A. H. Dayem and J. J. Wiegand, Phys. Rev. 155, 419 
(1967). 

G. M. Ehashberg, JETP Lett. 11, 114 (1970). 

V. T. Petrashov, V. N. Antonov, P. Delsing, and R. Glae- 
son, Phys. Rev. Lett. 70, 347 (1993). 

*° R. Shaikhaidarov, A. F. Volkov, H. Takayanagi, V. T. Pe- 
trashov, and P. Delsing, Phys. Rev. B 62, R14649 (2000). 

''^ C. Shelly, T. Higgs, C. Checkley, R. Shaikhaidarov, and 
V. T. Petrashov, in Condensed Matter and Materials 
Physics CMMP2011 Abstract Book (Manchester, 2011), p. 
168. 

M. H. S. Amin, A. N. Omelyanchouk, and A. M. Zagoskin, 
J. Low Temp. Phys. 27, 616 (2001). 



'^^ M. H. S. Amin, A. N. Omelyanchouk, A. Blais, A. Maassen 
van den Brink, G. Rose, T. Duty, and A. M. Zagoskin, 
Physica G 368, 310 (2002). 

M. Alidoust, G. Sewell, and J. Linder, Phys. Rev. B 85, 
144520 (2012). 

D. A. Ivanov and Y. V. Fominov, Phys. Rev. B 73, 214524 

(2006) . 

V. Braude and Y. V. Nazarov, Phys. Rev. Lett. 98, 077003 

(2007) . 

J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 

W. Belzig, F. K. Wilhelm, C. Bruder, G. Schon, and A. D. 

Zaikin, Superlattices and Microstructures 25, 1251 (1999). 

N. B. Kopnin, Theory of Nonequilibrium Superconductivity 

(Oxford University Press, Oxford, 2001). 

A. F. Volkov, A. V. Zaitsev, and T. M. Klapwijk, Physica 

C: Superconductivity 210, 21 (1993). 

Y. V. Nazarov, Phys. Rev. Lett. 73, 1420 (1994). 

A. V. Zaitsev, A. F. Volkov, S. W. D. Bailey, and C. J. 

Lambert, Phys. Rev. B 60, 3559 (1999). 
^'■^ E. Kandelaki, A. F. Volkov, K. B. Efetov, and V. T. Pe- 
trashov, Phys. Rev. B 82, 024502 (2010). 
^'^ K. D. Usadel, Phys. Rev. Lett. 25, 507 (1970). 

M. Y. Kupriyanov and V. F. Lukichev, Sov. Phys. JETP 

67, 1163 (1988). 

W. L. McMillan, Phys. Rev. 175, 537 (1968). 
^"^ L. G. Aslamazov and A. F. Volkov, in Nonequilibrium 
Superconductivity, edited by D. N. Langenberg and A. I. 
Larkin (Elsevier, Amsterdam, 1986). 

E. V. Bezuglyi, V. S. Shumeiko, and G. Wendin, Phys. 
Rev. B 68, 134506 (2003). 

I. O. Kulik and I. K. Janson, The Josephson Effect m 
Superconductive Tunneling Structures (Israel Program for 
Scientific Publications, Jerusalem, 1972). 

'''^ K. K. Likharev, Dynamics of Josephson Junctions and Cir- 
cuits (Gordon and Breach, New York, 1986). 

®^ A. Barone and G. Paterno, Physics and Applications of the 
Josephson Effect (John Wiley & Sons, New York, 1982). 
A. Y. Kitaev, Phys. Usp. 44, 131 (2001). 

F. Dominguez, F. Hassler, and G. Platero, e-print 
arXiv: 1202.0642 (2012). 

N. Didier, M. Gibertini, A. G. Moghaddam, J. Knig, and 
R. Fazio, e-print arXiv: 1202.6357 (2012). 
®^ P. Kotetes, A. Shnirman, and G. Schon, e-print 
arXiv:1207.2691 (2012). 

T. Yokoyama, Y. Tanaka, and A. A. Golubov, Phys. Rev. 
B 75, 134510 (2007). 

F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. 
B 68, 064513 (2003). 

A. F. Volkov, Y. V. Fominov, and K. B. Efetov, Phys. Rev. 
B 72, 184504 (2005). 

A. A. Golubov, M. Y. Kupriyanov, and Y. V. Fominov, 
JETP Letters 75, 588 (2002). 
^° V. N. Krivoruchko and E. A. Koshina, Phys. Rev. B 66, 
014521 (2002). 

T. Yokoyama, Y. Tanaka, and A. A. Golubov, Phys. Rev. 
B 72, 052512 (2005). 
'^^ T. Yokoyama, Y. Tanaka, and A. A. Golubov, Phys. Rev. 
B 73, 094501 (2006). 



